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We compare current and forecasted constraints on dynamical dark energy models from Type la 
supernovae and the cosmic microwave background using figures of merit based on the volume of 
the allowed dark energy parameter space. For a two-parameter dark energy equation of state that 
varies linearly with the scale factor, and assuming a flat universe, the area of the error ellipse can be 
reduced by a factor of ~ 10 relative to current constraints by future space-based supernova data and 
CMB measurements from the Planck satellite. If the dark energy equation of state is described by 
a more general basis of principal components, the expected improvement in volume-based figures of 
merit is much greater. While the forecasted precision for any single parameter is only a factor of 2-5 
smaller than current uncertainties, the constraints on dark energy models bounded by — 1 < w < 1 
improve for approximately 6 independent dark energy parameters resulting in a reduction of the total 
allowed volume of principal component parameter space by a factor of ~ 100. Typical quintessence 
models can be adequately described by just 2-3 of these parameters even given the precision of 
future data, leading to a more modest but still significant improvement. In addition to advances in 
supernova and CMB data, percent-level measurement of absolute distance and/or the expansion rate 
is required to ensure that dark energy constraints remain robust to variations in spatial curvature. 



I. INTRODUCTION 

In the absence of physically compelling models for dark 
energy, figures of merit (FoMs) are a useful tool that en- 
capsulate the constraining power of cosmological data. 
FoMs combine various constraints on the expansion his- 
tory of the universe (e.g. the distance-redshift relation) 
in a single number, or at most a handful of numbers, 
which serve as simple and quantifiable metrics by which 
to evaluate current and proposed experiments [ll-Q ■ The 
simplest schemes adopt fixed functional forms for the evo- 
lution of the dark energy equation of state and define the 
FoM as the inverse of the allowed parameter volume (if. 
One widely used version is the two dimensional WQ-w a 
parametrization [l| but other higher dimensional versions 
have also been considered 0-0]. 

Any simple parametrization of the expansion history 
risks biasing the FoM in favor of or against certain types 
of data by choosing a fixed functional form Q ■ To avoid 
this problem, one can use more complicated schemes that 
parametrize the whole functional freedom in the dark en- 
ergy equation of state evolution and separate the expan- 
sion history and growth of structure information. For 
example, uncorrelatcd modes of piecewise-constant dis- 
cretizations of the equation of state that are local in 
redshift I7H15H or constructed from principal components 



(PCs) [fl. [l6H2l| have been employed to characterize both 



current and future data. In particular, the inverse param- 
eter volume of the PC amplitudes, defined separately for 
each experiment, has been advocated as a FoM 0. 

Besides avoiding biasing results towards a particular 
functional form, a more model-independent FoM has the 
advantage of being able to identify improvements in dark 
energy constraints that might be missed by FoMs with 



fewer parameters. On the other hand, not all improve- 
ments in a multidimensional PC FoM reflect improve- 
ments in constraining the space of known dark energy 
models [IH, HU- For example, stronger constraints may 
exclude regions of the parameter space not occupied by 
typical models. 

In this paper, we study the FoMs defined both with the 
commonly-used wo-w a parametrization and with PCs 
based on forecasts for Type la supernova (SN) and cosmic 
microwave background (CMB) data. Previous studies of 
PC-based FoMs have generally relied on the Fisher ma- 
trix approximation, and the implementation and utility 
of these FoMs for real data has not been addressed. Here 
we define straightforward generalizations of PC FoMs 
and apply them to both actual data from current mea- 
surements and forecasts for future data. While we still 
construct the PCs using the Fisher matrix approach, we 
then treat the PC amplitudes as free parameters and 
constrain them using full Markov Chain Monte Carlo 
(MCMC) likelihood explorations with current or future 
data. This methodology follows that employed in previ- 
ous papers where we studied generalized predictions of 
classes of dark energy models based on forecasts [l7| and 
current data [l8[. The new element here is the applica- 
tion of these methods to the study of dark energy FoMs. 

The standard approach for constructing the principal 
component-based FoM is to use the PCs specific to the 
experiment and the cosmological probc(s) considered; sec 
e.g. Ref. Q. However, this approach makes it difficult 
to directly compare PC-based FoMs for different exper- 
iments and to assess improvements in specific regions of 
parameter space. To facilitate such direct comparison 
between FoMs for current and future data, we choose in- 
stead to fix the set of PCs based on a specific projection 
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FIG. 1: Redshift distributions of Type la supernovae in 
the Union compilation (solid blue) and assumed for SNAP 
(dashed red), including the anticipated low-redshift sample 
from other surveys at z < 0.1. 



0.15 - 



0.14 



A 

B 



0.13 



0.12 




13500 



14000 14500 
D(z.) [Mpc] 



of future data rather than computing separate PCs for 
future and current data. 

This paper is organized as follows. In Sec. [TTJ we de- 
scribe the current and forecasted SN and CMB data sets 
as well as additional priors from baryon acoustic oscilla- 
tions (BAO) and measurements of the Hubble constant. 
We also briefly review the MCMC analysis methods used 
to infer dark energy constraints. In Sec. IIII1 we com- 
pare FoMs from the current and future SN and CMB 
constraints for both the commonly-used wo~w a model 
(Sec. IIII A|) and a more general PC-based parametriza- 
tion (Sec. IIII 5|) . We summarize and discuss these results 
in Sec. [iVl 



II. METHODOLOGY 

In this section, we review the current data sets and 
assumptions about future experiments for forecasts that 
we use in this paper. We refer the reader to Refs. [13] 
(hereafter MHH1) and [H| (MHH2) for more details con- 
cerning the forecasts and current data sets, respectively. 
All forecasts in this paper assume that the data originates 
from a flat cosmological constant (w = —1) model with 
present matter fraction Q m = 0.24 and Hubble constant 
H a = 73 km s" 1 Mpc" 1 . 



A. Current SN and CMB Data 

The Type la SN sample we use is the Union compi- 
lation [2J|. These SN observations measure relative dis- 
tances, D(zi) I 'D{z2) 1 over a range of redshifts spanning 
0.015 < z < 1.551, with most SNe at z < 1 (see Fig.p. 
We include SN constraints using the likelihood code for 
the Union data sets (2f|, which includes estimated sys- 
tematic errors in the covariance matrix [24| . 

For the current CMB data, we use the 5-year data 
release from the WMAP satellite [26l - [28j employing the 
likelihood code available at the LAMBDA web site 12911 . 



FIG. 2: Approximate constraints on D(z*) and Qmh 2 from 
WMAP 5-year data (outer contours, blue shading) and fore- 
casts for Planck (inner contours, red shading), showing 68% 
CL (light shading) and 95% CL (dark shading) contours. 



Unlike the CMB priors used for the forecasts below, the 
likelihood used here contains the full information from 
the CMB angular power spectra, except for the small 
effects of gravitational lensing of the CMB that add lit- 
tle to current dark energy constraints from WMAP. We 
compute the CMB angular power spectra using the code 
CAMB [3(| H3 modified with the parametrized post- 
Friedmann (PPF) dark energy module [HI, [33| to include 
models with general dark energy equation of state evolu- 
tion where w(z) may cross w = — 1. 



B. SN and CMB Forecasts 

For our SN forecasts, we take the expected redshift dis- 
tribution for the Super Nova/ Acceleration Probe (SNAP) 
[llii] plus a low-z sample of 300 SNe at 0.03 < z < 0.1. 
The SNAP magnitude errors include both statistical and 
systematic components: 
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(1) 



where N a , shown in Fig. [TJ is the number of SNe in 
each bin of width Az (Az = 0.1 except for the sta- 
tistical uncertainties in the low-z SN bin, for which 
Az = 0.1 - zf n f n = 0.07), and Az sub is the width of 
the sub-bins used to smooth the distribution of SNe in 
redshift. We use 500 sub-bins up to z max = 1.7. The 
second term on the right hand side of Eq. (JlJ models a 
systematic floor that increases linearly with z up to a 
maximum at z max of 0.02 mag per Az = 0.1 bin [36| . 
For the Planck CMB forecasts, we use a 2 x 2 covariance 
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matrix corresponding to the inner error ellipses in Fig. [2] 

C CMB = 



(0.0018) 2 -(0.0014) 2 
-(0.0014) 2 (0.0011) 2 



with parameters 

g CMB = m^^/Mpc), n m h 2 } 



(3) 



Here D(z*) is the comoving angular diameter distance to 
recombination. 

In the Planck forecasts we ignore additional CMB in- 
formation about dark energy such as the ISW effect and 
gravitational lensing of the CMB. Planck data are ex- 
pected to obtain limits on the fraction of early dark 
energy of ct(Qde(z*)) ~ 0.004 [33, [38|; however, since 
these limits may depend on the modeling of early dark 
energy, we include in our forecasts only a conserva- 
tive prior of er(fiDE(z*)) = 0.025 (with fiducial value 
^de(z*) ~ 10~ 9 ) which approximates the current con- 
straint from WMAP [H. 



the simple WQ-w a models and general PC parametriza- 
tion. From the likelihood £(x|0) of the data x given 
each proposed parameter set 8, Bayes' Theorem tells us 
the posterior probability distribution of the parameter 
set given the data 



V(8\x) = 



£(x|0)7?(g) 
JdB C(x\0)V(oy 



(4) 



where V(8) is the prior probability density. The MCMC 
algorithm generates random draws from the posterior 
distribution. We test convergence of the samples to 
a stationary distribution that approximates 7^(0 Ix) by 
applying a conservative Gelman-Rubin criterion (47| of 
R — 1 < 0.01 across a minimum of four chains for each 
model class. We use the code CosmoMC [H, |49[ for the 
analysis of current data and an independent MCMC code 
for forecasts. 



III. FIGURE OF MERIT COMPARISONS 



C. Additional Priors 

For both the current and forecasted constraints we add 
additional priors from recent measurements of baryon 
acoustic oscillations and the Hubble constant. The BAO 
constraint we use is based on the measurement of the 
correlation function of SDSS Luminous Red Galaxies 
(LRGs) Jijj, which determines the distance and ex- 
pansion rate at Zbao ~ 0.35 through the combination 
D v (z) = [zD 2 (z)/H(z)} 1 / 3 . We implement this con- 
straint by taking the volume average of this quantity, 
(Dv), over the LRG redshifts, 0.16 < z < 0.47, and com- 
paring with the value of A = (Dy)\fVL m h 2 / zbao given in 
Ref. |4(|, A = 0.472±0.017 (taking the scalar spectral tilt 
to be n s = 0.96). Using more recent BAO constraints, 
e.g. from Ref. [4l|, has only a small effect on the current 
constraints with SN and CMB data [18[ . 

We include the recent Hubble constant (Ho) constraint 
from the SHOES team |42j, based on SN distances at 
0.023 < z < 0.1 that are linked to a maser-determined 
absolute distance using Cepheids observed in both the 
maser galaxy and nearby galaxies hosting Type la SNc. 
The SHOES measurement determines the absolute dis- 
tance to a mean SN redshift of = 0.04, which ef- 
fectively corresponds to a constraint on Hq for models 
with relatively smooth dark energy evolution in the re- 
cent past (cf. [43j]) such that lim^o D(z) = cz/Hq. Wc 
implement this constraint as a measurement of D(zh) = 
C2fc/(74.2±3.6 km s" 1 Mpc" 1 ). 



D. MCMC Methodology 

Given the current or forecasted data, we use MCMC 
likelihood analysis (e.g. see [44T4461 ] ) to determine dark en- 
ergy parameter constraints and figures of merit for both 



A. wq— w a Figure of Merit 

We first consider the t wo-p arameter model for the dark 
energy equation of state [50, [H| 



w(z) = w + ■ 



1 



(5) 



The FoM for this model defined by the Dark Energy Task 
Force [3,5 is the inverse of the area of the 95% CL region 
Ag 5 in the Wo-w a plane. Figure [3] shows these regions 
for the current and forecasted data with and without 
marginalization of spatial curvature J7k- 

For a Gaussian error distribution, Ag 5 is proportional 
to the square root of the determinant of the 2D covariance 
matrix C for wq and w a . Since the constant of propor- 
tionality used in practice for the FoM varies widely in the 
literature (e.g., see [521]). we simply define 
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(detC)- 1 / 2 

6.177T 
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(6) 



The approximate equality between the two lines in 
Eq. (jnj) becomes exact for a Gaussian posterior distribu- 
tion. Although the posterior in wq and w a is not perfectly 
Gaussian, the FoM computed using det C in Eq. ([6]) re- 
mains a good approximation to the area-based FoM. The 
difference between the two is ~ 10% in the worst case 
(current data with J7k ^ 0) and ^ 2% in the other cases. 
Values of the det C version are given m 

Table |U The FoMs for current data are consistent with 
those found in previous studies of WQ-w a constraints from 
similar data sets (e.g., [ol [53Tj55l |) . 

While analyses of WQ-w a models typically allow w(z) 
to cross —1, it is useful to also consider a more re- 
stricted class of models that satisfy the quintessence 
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FIG. 3: Constraints on wq and w a , assuming a flat uni- 
verse (left panel) or marginalizing over spatial curvature (right 
panel). 68% CL (light shading) and 95% CL (dark shading) 
regions are shown for both current Union+WMAP constraints 
(outer contours, blue shading) and SNAP+Planck forecasts 
(inner contours, red shading). Dashed lines mark the bound- 
ary of the quintessence prior, —1 < w(z) < 1. 



TABLE I: Figures of merit for wo~w a models and 10-PC 
quintessence models without early dark energy. 
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bound — 1 < w < 1 in light of the PC description be- 
low. Imposing this quintessence prior greatly restricts 
the allowed parameter space as shown in Fig. [3] Assum- 
ing a flat (nonflat) universe, the det C FoM is a factor 
of ~ 8 (6.5) larger than without the quintessence prior 
for current data, and a factor of ~ 13 (9) larger for fore- 
casts. Thus the addition of the quintessence prior in- 
creases the ratio of future to current FoM values by about 
60% (35%) relative to the ratios in Table |U Note that 
the effect of the quintessence prior on the FoM for fore- 
casts depends on the choice of the true model from which 
the data is drawn; for example, had we chosen fiducial 
wq and w a values that lie as far within the quintessence 
prior boundaries as the current data allow, the forecasted 
area allowed within the priors would be greater and the 
improvement in the FoM relative to current constraints 
would be smaller. 

Even in this case where the posterior distribution with 
the quintessence prior is far from Gaussian, det C still ap- 
proximates the area-based FoM through Eq. © reason- 
ably well. Although values of the det C FoM are smaller 
by 35-40% than the area FoM, the ratio of future to cur- 
rent FoM values is nearly unchanged. We shall see that 
this type of agreement for quintessence models carries 



over to to the more general PC-based FoM in the next 
section. 



B. Principal Component Figure of Merit 

We generalize the dark energy parametrization to allow 
arbitrary variations of the equation of state at z < z max 
with a basis of principal components (PCs) @, [U|. De- 
tails of the PC construction can be found in MHH1; here 
we highlight the points of special relevance for the FoM. 

We construct the PCs based on the SN and CMB fore- 
casts from Sec. Ill Bl The PCs are a set of orthogonal 
functions ordered by the precision with which they can be 
measured by the future SN and CMB data. Specifically, 
the principal component functions e%{zj) are eigenvectors 
of the SNAP+Planck covariance matrix in the space of 
piecewise constant values of the equation of state in red- 
shift bins, w(zj). The principal components form a basis 
in which an arbitrary function w(zj) may be expressed 
as 



w(Zj) = -1+ Yl 



(7) 



where at are the PC amplitudes, N Zj pc = 1 + ;z max /A,j 
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FIG. 4: Top panel: The first 10 PCs of w(z) (increasing vari- 
ance from bottom to top), with 500 redshift bins between 
z = and 2 max = 1.7. The PCs are offset vertically from 
each other for clarity. Bottom panel: an example of w(z) for 
a quintessence model (red dashed curve; see Eq. (|12[) ) and its 
representation using the first 10 PCs in Eq. (J7J (blue solid 
curve). Dotted lines show the w = — 1 zero point for each 
component (top) and for the example model (bottom). 
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FIG. 5: Marginalized ID posterior distributions for the first 6 PCs of fiat (solid blue curves) and nonflat (dashed red curves) 
quintessence models without early dark energy. Marginalizing over the early dark energy parameter Woo has little effect on the 
distributions. Top row: current Union+WMAP data; bottom row: forecasts for SNAP+Planck assuming a realization of the 
data with o»j = 0. Plot boundaries that cut off the distributions at nonzero probability correspond to the top hat priors on on 
for quintessence models. Higher variance PCs are mainly limited by the quintessence priors for both current constraints and 
forecasts (see Table HT]). Vertical dotted lines show the predictions of an example quintessence model from Eq. (|12l) (see Fig. [3]). 



TABLE II: Top hat prior rms Aoii/\/l2 and mean and rms of 
on from current data and forecasts, assuming flat quintessence 
models without early dark energy. 
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is the number of redshift bins of width Az, and Zj = 
(j — l)Az. We choose the maximum redshift for vari- 
ations in w(z) to be z max = 1-7, matching the largest 
redshift for the SNAP supernova data. The impact of 
dark energy evolution at higher rcdshifts is expected to be 
small, but perhaps no n- negligible; to account for this pos- 
sibility, we parametrize the equation of state at z > z max 
by a constant, Wqo. Likewise, we consider models with 
spatial curvature, f^K 7^ 0. Note that and fix are 
not allowed to vary from their fiducial values of —1 and 
0, respectively, in the Fisher analysis used to construct 
the PCs. 



Since the highest-variance PCs correspond to modes 
of w(z) to which even future data are insensitive, we 
truncate the sum in Eq. ([7]) by replacing iV Zj pc with 
N max < N Z , PC - As shown in MHH1, the 10 lowest- 
variance PCs (AT max = 10) form a basis which, for the 
classes of models we consider here, is sufficiently complete 
for the forecasts. Therefore, 10 PCs more than suffice for 
the current data as well. This set of basis functions is 
displayed in Fig. 01 

We impose priors on the PC amplitudes corresponding 
to the range of w(z) allowed for scalar field quintessence, 
— 1 < w < 1 following MHH1. Our conservative im- 
plementation excludes only parameter values that must 
violate these bounds even when considering possible com- 
pensation from the omitted higher-variance PCs (e.g. see 
the lower panel of Fig. 2] at z < 0.1). This approach 
yields top hat priors of width 
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(8) 



which follows from Eq. (A10) in MHH1. 

In analogy to the Wo—w a FoM in the previous section, 
we base the FoM for dark energy PCs on the determinant 
of the covariance matrix of a% for the n lowest- variance 
PCs, detC„. Even without informative data, this deter- 
minant is finite: 



det ci prior > = J 



Ac 



(9) 



where the individual factors of Aai/VT2 (listed in Ta- 
ble [n]) are the rms values of the corresponding top hat 
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( PC) 

FIG. 6: Top panel: PC figures of merit FoM„ with fore- 
casted uncertainties for SNAP+Planck and with measured 
uncertainties for Union+WMAP, normalized as in Eq. (|10[) 
to account for the quintessence prior. Bottom panel: Ra- 

f PC) 

tios of FoM„ forecasts to current values. In both panels, 
point types indicate different quintessence model classes: flat 
(solid points) or non-flat (open points), either with (squares) 
or without (circles) early dark energy. 



priors. These priors impact the FoM in a manner similar 
to those imposed in Ref. Q, where a Gaussian prior is 
adopted that requires the rms variation of 1 + Wi for each 
PC mode, averaged over scale factor, to be no more than 
unity. However, our priors are slightly stronger since we 
use the quintessence bounds to impose a top hat prior on 
1 + w(z) at all redshifts. In order that the FoM values do 
not reflect information that comes exclusively from the 
quintessence prior, we follow the convention in Ref. Q 
and rcnormalize the det C statistic to obtain 



FoMi PC ) EE 



det C„ 



detC 



(prior) 



-1/2 



(10) 



Table UU shows the mean and rms of each of the 10 
PCs from MCMC likelihood analysis using the current 
data and forecasts. Note that for components 6-10 in 
the current data the rms is dominated by the prior, and 
likewise for 7-10 in the forecast. Figure [5] shows the one- 
dimensional posterior probability distributions of the first 
6 PCs where the information from the data resides. 

Figure [5] (upper panel) shows the FoM for the current 
data and forecasts. Note that for both data sets the FoM 
starts to saturate around the 6th PC as expected. The 
lower panel shows the ratio of FoMs of future and current 
data. Once again the saturation point is around the 6th 
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FIG. 7: Forecasts for Qk and the amplitude of the first PC, 
ai, showing 68% CL (light shading) and 95% CL (dark shad- 
ing) regions. The large contours (blue shading) include only 
SNAP and Planck data as well as the additional priors from 
Sec. Ill Cl The small contours (red shading) add to these data 
a 1% measurement of (Dv) averaged over the redshift bin 
0.8 < z < 1.2 as might be obtained from a future BAO exper- 
iment, reducing the degeneracy between curvature and dark 
energy. 



PC with the total level of improvement varying from a 
factor of ~ 60 to ~ 140 depending mainly on whether 
spatial curvature is included. 

This dependence of the FoM improvement on curvature 
is mainly due to a degeneracy between fix and the first 
PC which reduces the FoM when marginalizing over 51k, 
particularly for the forecasts. This degeneracy is largely 
a consequence of our choice not to marginalize over J7k 
when constructing the PCs, but the leading degeneracy 
between dark energy and curvature is included in the 
first PC. There is a related difference between the flat 
and nonflat cases in the WQ-w a contours in Fig. [3] and 
FoM^"™"' in Table |U 

Current BAO and Ho measurements (Sec. Ill CP con- 
strain S7k enough that current SN+CMB PC uncertain- 
tics arc affected little by curvature, but for the SN+CMB 
forecasts we need ~ 1% measurements of an absolute dis- 
tance scale to achieve a FoM improvement comparable 
to that in the flat case. Figure [7] shows an example of 
breaking the ilx^cti degeneracy using a 1% (Dy) mea- 
surement at z = 0.8-1.2 as might be achieved from a 
future BAO experiment. As shown in Fig. 7 of MHH1, 
the shift in absolute distance corresponding to the curva- 
ture degeneracy is largely independent of redshift and its 
elimination could also be achieved with comparable mea- 
surements at z = — 0.1 from improved Hubble constant 
probes or z = 3 from high-rcdshift BAO. 

The small gap in ratios in Fig. [S] between 
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models with and without early dark energy is driven by 
the current constraint on the 5th PC; as shown in Fig.[5j 
P(ab) is not centered on zero and the distribution is cut 
off by the prior. Including early dark energy shifts this 
distribution further outside the prior, reducing the PC 
volume (and increasing the FoM) for current data. Fore- 
casts for PCs are basically unaffected by early dark en- 
ergy (see Fig. [5]), so the overall effect of early dark en- 
ergy is a slightly smaller FoM ratio when the 5th PC is 
included. This particular feature is not generic in the 
sense that it would not necessarily show up for other 
choices of data sets or dark energy parametrizations, but 
similar effects could appear in other analyses where the 
dark energy priors play an important role. 

Since we are using the PCs constructed assuming the 
forecasted data rather than the current data, we could 
in principle have strong covariances between the PC am- 
plitudes that would be hidden in Fig. [5] Likewise, we 
could have substantial differences between FoM^ PC ' ) and 
definitions of the PC FoM involving ratios of <Ji 0,0]. 
In practice, however, the covariances between PC ampli- 
tudes remain small in all cases. Even for current data, the 
difference between (detC„) -1 / 2 and n"=i ^ s ~ ^0% 
for flat quintessence models and < 30% for nonflat mod- 
els. The effect of covariances on the forecast FoM values 
is even smaller. Thus the ID distributions in Fig. [5] ac- 
curately depict the current and future constraints on PC 
amplitudes. 

A separate question is whether the use of det C gives 
misleading results due to non-Gaussianity of the PC pos- 
terior distributions (for example, due to distributions be- 
ing cut off by the quintessence prior). To test the sig- 
nificance of such effects, we consider an alternate FoM 
analogous to the area-based FoM for wo~w a models, 




where a is the parameter subset consisting of the first 
n PC amplitudes. The motivation of this form is that 
the allowed volume in parameter space is proportional 
to the inverse of the average number density of MCMC 
samples. The number density is proportional to V, and 
the averaging over the posterior probability gives another 
factor of V . The normalization of Eq. (fTTj) is chosen so 
that for an nD Gaussian posterior with covariance C„, 
V- 1 = (detC„)" 1/2 . 

For the 4-5 lowest-variance PCs, we find good agree- 
ment between U„ -1 and (detCn)" 1 / 2 with differences of 
no more than ~ 30%. Using the U„ _1 FoM increases the 
ratio of future to current FoM values by 10-15%. For 
n > 5, accurate computation of the integral in Eq. (fTTj) 
becomes more difficult due to the sparsity of MCMC sam- 
ples. 

Improvements in FoM^ PC ' do not necessarily represent 
significant improvements in the ability to limit the pa- 
rameter space of known dark energy models, especially 
for the higher PCs [22|, [23|. In order to address such is- 
sues, one can project the predictions of a model for w(z) 



onto the PC space and examine whether the predictions 
lie in the volume excluded by forecasted constraints. As 
an illustrative example, we consider a quintessence model 
with the potential 

V(<f>) = V + l -m 2 4> 2 . (12) 

This model provides examples in the thawing class (56j 
for Vo — > and can also mimic the low rcdshift behavior 
of Albrecht-Skordis models with oscillations around an 
offset minimum [57| . We consider an example where m = 
7 x 10" 33 eV and V /p cti t = 0.717 with fl DE = 0.733, 
h = 0.69 and Hubble-drag frozen initial conditions fa = 
0; the equation of state for this model is shown in the 
lower panel of Fig. @] These parameters are chosen to 
be allowed by the current data but testable with future 
data. 

Figure [5] compares the predictions for this model with 
the data constraints. Note that even for the future data 
only the first two PCs are stimulated at a level that 
the data can constrain. Thus for these types of mod- 
els FoM^ ' 1 is more representative of the parameter vol- 
ume improvements than higher-dimensional FoMs. Tak- 
ing into account the quintessence prior, the improvement 
from current to future data in this case is comparable to 
that implied by the WQ-w a FoM. 

More generally, studies have found that up to 3 PCs 
are useful for distinguishing amongst different commonly- 
used quintessence models with data sets comparable to 
our future forecasts 23] . On the other hand, these studies 
do not preclude the possibility that FoM improvements 
in the higher components can distinguish between other 
yet to be investigated classes of models. 

IV. DISCUSSION 

We have compared current and forecasted figures of 
merit (FoMs) for dark energy using both a simple Wo—w a 
description of the equation of state, and a more com- 
plicated but more complete principal component (PC) 
paramctrization. By consistently using a fixed set of PCs 
based on future data, and by generalizing the FoM defini- 
tion to include possible parameter covariance, we showed 
how PC FoMs can be applied to likelihood analyses of 
both existing data sets and forecasts. 

We have also shown that the covariance-based FoMs 
accurately represent relative changes in the parameter 
volume, even in the presence of strongly non-Gaussian 
posterior distributions such as those caused by impos- 
ing top hat priors. Traditional variance-based PC FoMs 
are consistent with those that include the full covariance 
for the cases we have tested, but we do not expect such 
agreement to hold in general. For example, if our PC 
eigenfunctions from SN and CMB data are applied to 
qualitatively different data, e.g. weak lensing or BAO, the 
FoMs that account for the parameter covariances should 
be employed. 
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For the WQ-w a FoM, future space-based supernova 
data and CMB measurements from the Planck satel- 
lite can improve on current measurements by a factor 
of ~ 10. For the PC FoM, the expected improvement 
is much greater still. While the forecasted precision for 
any single principal component is only a factor of 2-5 
smaller than current uncertainties, the constraints on 
general quintessence models bounded by — 1 < w < 1 
improve for approximately 6 components resulting in a 
reduction of the total allowed volume of dark energy pa- 
rameter space by a factor of ~ 60-140. The FoM im- 
provement depends mainly on whether or not variations 
in spatial curvature are allowed, with the maximum ratio 
of 140 requiring cither fixing the curvature with a the- 
oretical prior (fix = 0) or combining future SN and CMB 
data with a measurement of absolute distance and/or the 
expansion rate with at least ~ 1% accuracy. 

Although improvements from future SN and CMB 
data extend to 6 components of w(z), many commonly- 



considered quintessence models are adequately described 
by the first three PCs for which the allowed volume only 
decreases by a factor of ~ 20-40 relative to current mea- 
surements. While figures of merit provide a convenient 
model-independent metric to evaluate experiments, they 
must be interpreted with caution when science goals in- 
clude discriminating amongst specific dark energy mod- 
els. 
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